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ABSTRACT 

We present near-infrared (2.5 — 5.0 /im) spectra of shocked H2 gas in the supernova remnant IC 
443, obtained with the sateUite AKARI. Three shocked clumps — known as B, C, and G — and one 
background region were observed, and only H2 emission lines were detected. Except the clump B, the 
extinction-corrected level population shows the ortho-to-para ratio of ^ 3.0. From the level popula- 
tion of the clumps C and G — both AKARPs only and the one extended with previous mid-infrared 
observations — we found that the v = levels are more populated than the f > levels at a fixed level 
energy, which cannot be reproduced by any combination of H2 gas in Local Thermodynamic Equilib- 
rium. The populations are described by the two-density power-law thermal admixture model, revised 
to include the collisions with H atoms. We attributed the lower (n(H2)=10^'®~'^'* cm~'^) and higher 
(n(H2)=10^'*~^'^ cm~^) density gases to the shocked H2 gas behind C-type and J-type shocks, respec- 
tively, based on several arguments including the obtained high H I abundance n(H I)/ri(H2)=0.01. 
Under the hierarchical picture of molecular clouds, the C-type and J-type shocks likely propagate 
into "clumps" and "clouds" (interclump media), respectively. The power-law index b of 1.6 and 3.5, 
mainly determined by the lower density gas, is attributed to the shock-velocity diversity, which may 
be a natural result during shock-cloud interactions. According to our results, H2 u = 1 ^ S(l) 
emission is mainly from J-shocks propagating into interclump media. The H2 emission was also de- 
tected at the background region, and this diffuse H2 emission may originate from collisional process 
in addition to the ultraviolet photon pumping. 

Subject headings: Infrared: ISM — ISM: individual (SNR IC 443) — (ISM:) supernova remnants — 
Shock waves 



1. INTRODUCTION 

IC 443 (G189.1-f3.0) is an extensively studied su- 
pernova remnant (SNR) and famous for its interac- 
tion with nearby molecular clouds. The interaction 
has been widely observed over div erse wavelength s: the 
7-ray from hardronic collisions (lEsposito et al] 
Albert et ani2007t [Acciari et al.|[2009HAbdo et"al1 



1996 



2010; 



Tavani et al.l !2O10lr _the X-ray absorption by fore- 
ground clo uds (iTroia e t al. 200 ^, t he infrared forbid- 
den lines (Burton ct al. 1988, 1990l; llnoue et all 119931: 
Richter et al. .1995a. : Ce sarsky et all 119991: iRho et al l 



ra- 



20011: IRosado et all [20071) . the enh anced CO line 
tio (|Seta et al.lll998|: IXu et al.ll201lD. the bro ad molec 



ular emission lines Jwhite et al.lll987t IWang fc Scovillc 
1992"; 'Dickman et al.' '1992': 'van Dishoe ck et al.l il993c 
SncU ct al. 2005; Zhang ct al. 2009), the bow-like 
featur e in position-velo city diagrams of molecular 
line (iTauber et al.l Il994f). and the OH maser line 



(IClaiissen et al.lll997i: IHoffman et al.l 120031 : iHewitt et al.l 
I2006D . Therefore, IC 443 is usually observed when study- 
ing the shock-cloud interaction. 

Its age ranges f rom ~ 3 — 4 kyr (jPetre et al.l 119881 : 
IWang et al][l99a iTroia et al.ll2008l) to ^ 20 - 30 kvr 
(|Chevalierill999l : l01bert et al.ll2001HGaensler et al.]l2006l: 
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iBvkov et al] 120081: iLee et al.l l2008aD . and its distance 
is thought to be about 1.5 kpc based on several ar- 
guments, such as the contact w ith Gem OB 1 as- 
sociation (IPoveda fc Woltieil Il968f). the empirica l S — 
D relatioinMilnel 119791: iCaswell fc Ler chd [T979l|. the 
remnant energy (Fesen fc KirshneH 119801 : iFesenl 
and the high-velocity absorption lines observed 
against background s tars (Welsh fc Sallmen .2003) . It 
extends - 45' (cf. IGaensler et al.l I2006D and over- 
laps in the sky with another rtiore e xtended SNR 
G189.6-h3.3 (lAsaoka fc Aschenbachlll99l . IC 443 con- 
sists of two half shells, and another large shell of the 
SNR G189.6-I-3.3 overlaps with the former two shells; 
these thr ee shells were named a s A, B, and C, respec- 
tively, bv lBraun fc Stroiiil (|1986( ). Th e overall pic t ure o f 
the remnant region is well outlined in lTroia et aH (|2006l ) 
andl Lee et"al] |2008a '). In the middle of the two shells A 
and B, there is a W-shaped r idge which show s strong H2 
emission lines (cf. Fig. ID and IRho et al.ll200il ): this ridge 
is thought to be a torus-type molecular clouds overrun 
by the SNR shock. 

Infrared H2 emission lines are useful to study shocked 
molecular clouds, since H2 is the most abundant molecule 
and its quantum levels cover a wide energy range enough 
to study shocked gas, whose temperature ranges from 
a few hundred to a few thousand kelvin. Toward 
the W-shaped H2 ridge, infrared spectra l observations 
have a lready been performed from ground (jRichter et al.l 
ll995 aD and space by Infrared Space Observatory (ISO, 
Ces arskv et all 119991 ) and Spitzer (Ncufcld ct al'. 2007). 
However, the H2 emission lines which falls within 2.5 — 5.0 
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fim have not been observed completely; the ground obser- 
vations missed several lines because of the atmospheric 
absorption, and this wavelength range is not covered by 
Spitzer spectroscopy and was simply not observed by 
Infrared Space Observatory (ISO). Besides, this wave- 
length range is worthwhile to observe because, for in- 
stance, we could obtain the population of high- J v ~ 
levels, which is usually assumed to follow v > lev- 
els (e.g. iRho et al.|[200lHGiannini et al.ll2006[) . but have 
not been thoroughly checked yet. 

Here we present the results of near-infrared spectral 
observations over 2.5 — 5.0 /im for the shocked H2 gas 
in the SNR IC 443. The observations were performed 
with the satellite AKARI (section[2]). We detected many 
H2 emission lines toward shocked molecular gas (section 
13. and found that the population of the shocked H2 
gas cannot be described by any combination of H2 gas in 
Local Thermodynamic Equilibrium (LTE), which have 
been usually used for (section [5^ and [53]) . I nstea d, we 
employed a non-LTE H2 gas model (section 13.41) , and 
interpreted the results in terms of a shock combination 
(C-type and J-type) together with the diversity in shock 
velocities (section l4.ip . The observed background emis- 
sion, attributed to diffuse warm H2 gas, was discussed as 
weh (section 1321) ■ 

2. OBSERVATIONS AND DATA REDUCTION 

The spectral observa tions were performe d with the In- 
fraRed Camera (IRC, ' Onak a et a P '2007D onboard the 
Japanese satellite AKARI (jMurakami et al.l [2007D . on 
2008-Sep-26th and 27th during the post-Helium phase. 
During this phase, only near-infrared observations were 
possible, since the cryogenic cooling with liquid Helium 
had been run out. We used the 5" x 48" slit and the 
grism, whose resolving power and wavelength coverage 
are AA ^ 0.03 /im and 2.5 — 5.0 /Ltm, respectively. The 
observation mode, called as Astronomical Observation 
Template (AOT), is IRCZ4 which is designed for gen- 
eral spectroscopic observations. It has an imaging obser- 
vation san dwiched by spectro scopic observations of four 
frames (cf. lOnaka et al.i |2009D. Comparing this refer ence 
image to the 2MASS catalog (jSkrutskie et al.ll2006D . we 
corrected the default astrometry of the slit, given from 
the satellite attitude information. 

We observed four regions, three of which are the 
shocked CO clumps and the rest is the background. Fig- 
ure [T] shows the four slit position s over the 2MASS Kg 
RGB image (jSkrutskie et al.|[2006( ). which displays a dif- 
fuse 'W' feature th at traces the Ha v = 1 ^0 S(l) 2.12 
/zm line emissions (IRho et al.ll200l . We name three on- 
source positions as 'B ', 'C, and 'G' after the names of the 
shocked CO clumps ()Denovei] [r979: Huang ct al. 1986), 
and the background as 'BG.' Table [1] summarizes our 
observations — the region name, RA-Dec position, obser- 
vation ID, and AOT. 

The data were reduced through t he official pipeline , 
supported by the AKARI team (cf. lOnaka et allboogf ). 
We used the new spectral response curve for the post- 
Helium phase datEQ which shows a degraded sensitivity, 
^70% of the Helium phase sensitivity. Columns of the 
two dimensional spectral images are occasionally satu- 
rated by the very bright sources in the imaging area of 

* http:/ /www. ir.isas.jaxa.jp/AKARI/Obscrvation/DataReduction, 



the detector, which cause the column pull-down effect. 
This effect was corrected by masking out the relevant 
columns. Hot pixels of the detector were also masked 
out. During the data reduction, no smoothing and tilt- 
correction were applied to the two dimensional spectral 
images. 

In order to extract spectra, we chose certain sections 
along the slit length, then averaged the pixel values 
within those sections. The extraction sections were care- 
fully chosen for the H2 emission lines to be nearly uniform 
within the sections. The sections extend 10 pixels for the 
clumps B and C, 14 pixels for the clump G, 25 pixels for 
BG, where one pixel corresponds to 1.46" (jOnaka et al.l 
[2003); for the data set 1420806-002 of BG, we only chose 
5 pixels as an extraction section to avoid abnormal sin- 
gle pixel peak. Figure [2] displays the extracted spectra 
at each region. For clarity, the error bars were omitted 
and the spectrum of BG was enlarged by a factor of 20. 
The statistical error bars can be seen in Figures [5][ni 

For the error estimation, we included the systematic er- 
ror caused by the calibration source type, in addition to 
the statistical error. Our target is a diffuse source, hence 
the flux calibration referred to standard point sources 
are not suitable for our source, since the aperture loss 
and the slit loss would vary with the source type; the 
official pipeline uses the calibration from point sources. 
We considered this type of systematic error and adopted 
10% of the signal intensity as the systematic error (pri- 
vate communication with the AKARI helpdesk). It was 
squarely summed to the statistical error of the line in- 
tensity, cr^( -I- cr^ys , after measuring the line intensities 
(cf. section [3]). 

3. ANALYSIS AND RESULTS 
3.1. Line Identification and Intensity Measurement 

Figure[2]shows many emission lines. Since our concerns 
are the H2 emission lines, we first compared the wave- 
length of the observed emission lines with those of H2. 
For easier identification of single and blended lines, we 
made template spectra of H2 gas in LTE at diverse tem- 
perature from 1000 K to 4000 K. This temperature range 
was adopted because the shocked H2 gas usually shows 
such a range of excitation temperatures at the upper lev- 
els of E(v, J) ~ 5x 10^-2. 5x 10^ K (e.g. IRosenthal et al.l 
l2000t IRho et al.l [200TI : iGiannini et alTl2006D . which in- 
cludes the upper levels of the H2 emission lines we de- 
tected. In this way, we identified all the detected emis- 
sion lines as H2 emission lines. 

Some lines are blended with nearby lines, hence we 
could not tag them as a single line. The uniquely iden- 
tified lines and the blended lines are indicated by '|' 
and respectively, in Figure [21 and their line iden- 
tifications are listed in Table [2] For a cross-check, this 
identification was compared with that of the shocked H2 
gas observed in the O rion Molecular Cloud-1 (OMC-1, 
IRosenthal et al.|[2000[ ). and our identification turned out 
to be reliable. Some contribution from Br/3 2.63 /im can 
be blended with the 2.63 /im blended line (cf. Fig. [2]); 
however, we think it is unlikely since Bra 4.05 /im, 
which should be stronger than Br/3, was not detected; 
the spontaneous transition probabilities of Bra and Br/3 
are A54 = 2.7 x lO*" s"^ and Aq^ = 7.7 x 10^ s~\ respec- 
^fj^JpRlsfeeilSafjaifefeS^een at the edge of the band (< 2.6 
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and > 4.9 fim) were ignored, since they are likely to be 
inadequate to analyze. The 2.56 /im and 4.95 /xm fea- 
tures seem to be the u = 1 — Q(9) blended with nearby 
lines and the u = 1 — 1 S(9), respectively. 

The line intensities were measured by fitting their 
line profiles with a continuum plus Gaussian whose fuU- 
width-at-half-maximum (FWHM) is fixed to the spectral 
resolution of the IRC (AA — 0.03 fj,m). Adjacent lines 
whose profiles are overlapped with each other were fitted 
simultaneously. As a baseline continuum for the fitting, 
we used a median-smoothed spectrum of each region; the 
kernel width ranging from 0.20 to 0.54 /im was carefully 
chosen to be wide enough to erase out the line emis- 
sion features. The feature near the edge of the spectrum 
(e.g. 2.56 /im feature in Figure |3]) remains unchanged af- 
ter the median filtering, since the filter only works on 
those pixels that are away from the edge by more than a 
half of the filter width. The fitted profiles are displayed 
in Figures [SHni and their measured intensities are listed 
in Table [2l The intensities of blended lines are shown 
with the symbol '<', since we cannot determine the indi- 
vidual contribution of each line blended. The intensities 
of weak lines, whose signal-to-noise ratio are lower than 
3.0, are expressed as 90% confidence upper limits. 

3.2. Reddening Correction and H2 Level Population 

Over 2.5 — 5.0 /im wavelengths, the extinction opti- 
cal depth becomes greater than one with the hydrogen 
nuclei column dens ity, N^=N{}1 l)+2N(H2), of > lO^^ 
cm' -2 (IDrainell2003l ). Since the observed regions are per- 
vaded with dense molecular gas, the measured intensity 
should be corrected for the reddening by the interven- 
ing interstellar dust. We e xploited the extinction curve 
of "Milky W ay, Ry = 3.1" (IWeingartner fc Drainel[200l 
iDraind [20031 ) ■ and adopted proper hydrogen nuclei col- 
umn densities for ea ch region. Av= 1 3.5 w as adopted 
for the clumps C as 'Neufeld & Yuad (|2008D did based 
on the results of Richtcr et al. (1995a)- We adopted the 
same Ay for the clump B, since its e xtinction is known t o 
be similar with that of the clump C (jBurton et al.lll988l ). 
No extinction measurement exist toward EG; thus, we as- 
sumed it to be the Ay of the clump C, the nearest clump 
from BG. ^^=10.8 was adopted for the clump G; this 
was inferred from A2.12 — 1.3, obtained bv lRichter et al.l 
(|T995a' ). employing the "Milky Way, i?y=3.1" extinc- 
tion curve. The corresponding iVn was calculated with 
the equation Av= A^h/(1-87 x lO^^ cm'^) for Rv=SA 
(|Bohlin et al.l[T97l . 

The H2 level population was derived from these red- 
dening corrected intensities, assuming that the II2 emis- 
sion lines are optically thin. Since the infrared H2 emis- 
sion is emanating from electric quadrupole transition, it 
is optically thin under a typical interstellar medium con- 
dition; for instance, the pure- rotational S(0) and 8(1) 
lines become optically thick at the line center when 
A^(Il2) > lO^'' cm~2, adopting the line width of 10 km 
s~^. We derived the reddening corrected level population 
of II2 gas from the following equation. 



Nrciv,J) 



47rA Irc{v, J v' , J') 

he A{v, J v', J') 



(1) 



where Irc{v,J — J> v',J') and A{v,J — ?• v',J') are the 
reddening corrected line intensity and the Einstein-A ra- 



diative transition probability of the transition from level 
(u, J) to {v',J'), respectively. The molecular data for 
H2 were obtained from the database provi ded by a sim- 
ulation code, CLOUDY (version C08.00; iFerland et all 
[T99i l. The results are listed in Table [3] and their popu- 
lation diagrams are displayed in Figure [71 gj is a weight 
factor which corresponds to (2 J -I- 1) and 3(2 J -I- 1) for 
para (even J) and ortho (odd J) states, respectively; the 
population of LTE gas is appeared as a straight line in 
this diagram. 

In Figure [II the clumps B, C, and G show the popula- 
tion of w = 0, 1, 2, while BG shows that of {v, J)=(0,11) 
only. In the clumps B, C, and G, the f = population 
shows a similar shape and little zigzag pattern; when 
the ortho-to-para ratio is approa ching to 3. , the zigzag 
pattern disappears (cf. Ncufcld eral1l2006L [20071 ). The 
V = 1 population, however, shows a similar shape and lit- 
tle zigzag pattern only in the clumps C and G; the clump 
B shows an evident zigzag pattern over the population 
of (u, J)=(l,l), (1,2), and (1,3). Here we note that the 
u = and v — I levels follow different branches. This 
becomes clearer when plotted with the lower- J v — 
population obtained from mid-infrared observations (see 
the following section). BG show a much smaller but ev- 
ident population of the {v, J)=(0,11) level. It is about a 
factor of 12 — 38 smaller than the clumps B, C, and G 
(cf. Table [3]). 

3.3. Comparison of H2 Populations with Previous 
Observations: the Clumps C and G 

As far as we k now, t he UKIRT CGS4 observation 
of iRichter et al.l (|1995aD is the only published near- 
infrared spectroscopic observation that targeted on the 
shocked II2 gas in IC 443; they observed the clumps 
C and G. Besides, toward these two clumps, there are 
published results of mid-infrared spectroscopic observa- 
tions for the shocke d II2 gas, performed with Infrared 
Space Observatory (jCesarskv et all |1999[ ) and Spitzer 
(jNeufeld et all |2007| ). For the clump B, mid- infrared 
spectros copic observations were al so performed with 
Spitzer (jNoriega-Crespo et al.ll2009[ l. however, they are 
under analysis and only covered three emission lines, 
w = - S(0), S(l), and S(2). Therefore, we concen- 
trated on the data of the clumps C and G, and compared 
our AKARI results with the ones from previous studies. 

Before the comparison of the results, we compared 
the spectra extraction areas of AKARPs and other's, 
to check whether we are comparing the same radiation 
sources. The yellow circle and boxes in Figure [S] indi- 
cate the extraction areas of the clumps C and G. The 
details for the extraction area are given in the figure 
caption. For the clump C, the centers of the extraction 
areas falls within ~ 20". Especially, that of Spitzer data 
(jNeufeld et al.l[2007f l — the center of a Gaussian taper — 
well falls into the extraction area of the AKARI data; 
hence, it is likely that the H2 emissions, observed by 
AKARI and Spitzer, represent almost the same gas, un- 
less the average physical property drastically varies over 
a few arcsec scale (< 0.07 pc). For the clump G, the 
centers of the extraction areas falls within about 30". 
In thi s case, the extr action areas of AKARI and ISO 
(Cesa rskv et al.|[T999D data overlap about 40 % with a 
~ 10" (~ 0.07 pc) separation. Thus, as in the clump C, 
the H2 emissions, observed by AKARI and ISO, likely 
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represent almost the same gas, unless the average physi- 
cal property drastically varies over a ~ 0.07 pc scale. 

Prior to the comparison of level population from differ- 
ent observations, we here note that the same extinction- 
correction curve and Ay values, used in section 13.21 
were used for the consistent comparison. We first com- 
pared oui^j^suhs_mtli_th observa- 
tions (iCesarskv et al.lll999l:lNeufeld et al.ll2007[) . The top 
panels of Figure [9] shows the comparison. We adopted 
15% systematic calibration error, which is a do minant 
component, for both clumps, C and G (Ccsarsk v et al.l 
119991: iNeufeld et ani2007f ). We here note that, in both 
clumps, V = levels are more populated than u = 1 
levels at a fixed level energy with a 3-cr significance at 
least. This means that the H2 level population cannot 
be reproduced by any combination of H2 gas in LTE, 
which is us ually adopted for the description of shocked 
H2 gas fe.g. lRho et aHIMlt iGiannini et al.l[2006 ). This 
invalidity was previously pointed out in the study on 
the shocked H2 g as of the supernova remnant HB 21 
(|Shinn et al.l[20ia) . 

Secondly, we compared our results with those of ground 
near-infrared observations by Richtcr ct al. (1995a). The 
extinction-corrected level populations were derived from 
their "averaged" intensities of "position 1" and "posi- 
tion 3," each for the clumps C and G, respectively. The 
comparisons are plotted in the middle panels of Figure 
[9l As the figures show, there is an inconsistency be- 
tween the same {v, J) levels of two observations. We 
attribute this inconsistency to the flux calibration differ- 
ence; the V — 0,1 levels show almost constant, vertical 
gaps in the population diagram, about a factor of 3 — 4 
in the column density. Our AKARI calibration is likely 
to be correct, because the AKARI v — level popu- 
lations are seamlessly merged with those obtained from 
previous mid-infrared observations (cf. the top panels 
of Fig. in]). However, we cannot rule out that the incon- 
sistency is caused by the difference of observed regions, 
since the average properties of the shocked H2 gas may 
change over ~ 30" ^ 0.2 pc scale (see Fig. [8|). 

Lastly, we note the importance of the space observa- 
tions, fully covering 2.5 — 5.0 /xm, for the study of shocked 
H2 gas. With ground observations, the gap between 
u = and 1 levels is hardly inspectable, since only a 
few V = levels can be observable at largely separated 
E(v, J) (see the grey plots in middle panels of Figure|9]). 
Hence, when using the results of ground observations, 
wc arc likely to think that a combination of H2 gas in 
LTE looks viable to reproduce the level populations up 
to E{v, J) ^ 25,000 K (see the bottom panels of Figure 
[9]). However, such a gap can be immediately inspectable 
from the space observations, fully covering 2.5 — 5.0 /im, 
since u = 0, 1 levels are covered continuously; the AKARI 
near-infrared observation is a good example. 

3.4. Power-law Thermal Admixture Model of II2 Gas 

As seen in the previous section, the combination of 
H2 gas in LTE cannot reproduce the observed level pop- 
ulation of the shocked H2 gas in the clump C and G, 
obtained from AKARI and previous mid-infrared ob- 
servations. Therefore, we applied the power-law ther- 
mal admixture model of H2 gas, which successfully re- 
produced the level populati on of shocked H2 gas before 
((Neufeld fc Yuan .2008. : .Shinn et al.ll200l iNeufcld et al.i 



200 1 iShinn et al.l [20I0I: iLee et al.l [20I0I: iTakami et ahl 
2010: 1 Yuan fc Neufeldl 120111) . The model configuration 



was the same with that use d in th e study of the super- 
nova remnant HB 21 (Shin n et al1 [2"009.^201 Ci) and th e 
young stellar object (YSO) L 1251A f|Lee et all IMol ) . 
except the inclusion of the H atom as an additional col- 
lider. Our AKARI spectra include several ro- vibrational 
transition lines (e.g. v = 1 ^ 0,2 1), which are sensi- 
tive to the collision with H atoms (jNeufeld fc Yuan|[2008l : 
IShinn et al.ll2009ll20Tol) . Thus, we updated the previous 
model to reflect the collisions with H atoms. 

The H2 column density was calculated from the follow- 
ing equation. 



dN = aT-^dT, 
e, a ^ 



(2) 



where, a = ''^"f^Z^Jr'^ 



T„ 



100 K, T„ 



: 4,000 K 



a and h are constants, and N(iJ2;T > 100 K) is a to- 
tal column density of molecular hydrogen warmer than 
100 K. At each temperature, the statistical equilibrium 
was assumed, and the coUisional partners were H2, He, 
and H. The coUisional dccxciation rates were obtained 
from iLe Bourlot et al.l ([1999 ) for H2 and He, and from 
iWrathmall et al.l()2007D for H. For th e H2 col lider, newer 
rates were calculated by ILee et ahl (l2008bD: ho wever, 
their results are similar with ILe Bourlot et aLl 's over 
100 — 6000 K and onl y include for those levels of J„p < 8. 
Hence, we kept using iLe Bourlot et al.l 's. The coUisional 
excitation rates were calculated from the detailed bal- 
ance relation. The H2 density, n(H2), and the relative 



n(H l) 
n(H2) 



were 



abundance of H atom to H2, Xh = log 

set as free parameters, while n(He) was assumed to be 
0.2xn(H2). The ortho-to-para ratio was set to 3.0, since 
the level populations show little zigzag pattern (cf. the 
top panels of Fig. 

First, we tried to fit with a single n(H2) value, but we 
failed. No single n(H2) model could successfully repro- 
duce the population of high energy-levels. This result 
reafhrms the previously-noted tendency that the single 
n(H2) model does not reproduce the level population over 
E{v,J) = - 25,000 K (cf. section SXl]) . Neither the 
addition of another b was successful. Instead, we added 
an additional n(H2) component, and then we obtained 
successful results. The fitting results are displayed in 
Figure [TOl For the clump G, the fewer data points from 
mid-infrared observations caused weaker constraints for 
the fitting parameters than the clump C; therefore, we 
fixed the Xh value as —1.7 in view of the fitting results 
for the clump C. Then, we scanned the space for fit- 
ting parameters and the results are displayed in Figure 
[TTl the step size for the scan was 0.1 for all parameters 
displayed in Figure ill! The fit parameters with their 
90% confidence intervals, together with the reduced chi- 
square values, are listed in Table S] 

As seen from Figure [TOl the reddening-corrected level 
populations are well reproduced by the power-law ther- 
mal admixture model with two different n(H2)s: one 
low n(H2) ~ lO^^ — 10^ cm~'^ and the other high n(H2) 
~ 10^ - 10^ cni-3 (cf. Table HI). The lower and higher 
density H2 gases mainly contribute to the lower and 
higher upper-energy levels, respectively. This indicates 
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that some model parameters may not be uniquely deter- 
mined, since only a small portion of the modeled popu- 
lation is constrained by the observed population. This 
point is more discussed in section l4.1.1l The two kinds of 
H2 gases share a common power-law index 6; it is 1.6 and 
3.5 for the clumps C and G, respectively. The column 
density ratios of the lower to higher density gas are not 
much different; they are 40 and 13 for the clumps C and 
G, respectively (see section lT.l.ll for more notes about the 
column densities, however). We obtained Xh = —1-7, 
which corresponds to ri(H I)/n(H2)=0.02; this value is 
similar with the one where the collision with H atoms 
starts to dominate the collision with H 2 for the rovibra- 
tional trans ition lines, a s ment i oned by [Neufeld &: YuanI 
(|2008D and iShinn et all (|2009l I2O1O0 . and smaller than 
those obtained for protos tellar outflows of LDN 1157, 
n(H I)/n(H2) = 0.1 - 0.3 (|Nisini et al.ll201fln . 

The x'^ contour for the column densities show a rough 
correlation for both clumps, C and G (cf. Fig. [TT|). This 
seems to be caused by the fact as follows. If both column 
densities increase or decrease together, then the can 
be decreased by adjusting the power-law index h shared 
by the lower and higher density H2 gases. However, if 
one column density increases and the other decreases, or 
vice versa, then the cannot be decreased by the same 
way as before, since the total H2 population changes its 
shape, which is related with the power-law index b. The 
variation of n(Il2) and iV(H2) are not effective as of b 
in adjusting the fitting, as the confidence intervals in- 
dicate (Table 2]). For the clump C, almost no correla- 
tion is shown between the densities and between b and 
Xh- For the clump G, the densities show a complex 
surface. It may be caused by the absence of the low- 
est level data, (u, J)=(0,2) and (0,3), which imposes a 
weaker constraint for the fitting. 

4. DISCUSSION 

4.1. Shocked H2 Gas Described by the Power-law 
Thermal Admixture Model 

4.1.1. Comparison Between Model Parameters from 
Previous and Our Studies 

The first attempt to describe the level population of 
shocked H2 g as with the power-l aw mixtu re of thermal 
H2 ga s was bv lBrand etaLllflQSl ( see also lBurton et al.l 
Il989| ). The y assumed that the H2 gas, shocked by J- 
type shocks (jPraine fc McKcc 1993), is in LTE and the 
postshock temperature profile is determined through the 
H2 radiative cooling, which is proportional to T^'^. In 
this way, they modeled the level population of thermally 
mixed H2 gas at the postshock region, which is equivalent 
to the power-law thermal admixture of H2 gas in LTE 
with the power-law index 6 = 4.7 (cf. section IX4|) . They 
applied this model to the observational data of OMC-1, 
and found that it is successful in describing the available 
data at that time. However, the omission of important 
coolants — like CO, OH, H2O, and grain — at dense en- 
vironments as OMC-1 and of the magnetic field make 
the model assumption doubtful (Chang fc M artin 199_li; 
IDraine fc McKea[T99l. 

Then. lOliva et aLf i 19901) first tried to use the method 
of power-law thermal mixing as a phenomenological de- 
scription tool for the H2 level population, with no back- 
ground physics; they simply mixed H2 gas in LTE ac- 



cording to the power- law distribution, dN{T) ~ T-^dT. 
In this approach, they found that the H2 population of 
the SNR RCW 103, obtained from near-infrared observa- 
tions, can be described with the models whose power-law 
indices are between 6 = 3.8 and b = 4 .7. This method 
was extended to the non-LTE case bv INeufeld fc YuanI 
()2008D . Applying to the Spitzer IRAC data, they found 
that the shocked H2 gas of the supernova remnant IC 443 
can be described with the power-law index & of ~ 3 — 6 
and the molecular hydrogen density of ri(H2) ~ 10^ — 10^ 
cm~^. Afterwards, the same meth od was a pplied to 
the s hocked H2 of the SNR HB 21 (IShinn et al. 200 " 



201c ) and of the outflows from YSOs ( Neufeld et al 



2009 : IShinn eTaTI IMl lLee et al.l [2010 



2oT0t INisini et al.l [20T(i[ lYuan fc Neufelc 



Takami et al.l 
I2011D. All of 



the non-LTE application res ults are summarize d in Table 
[5l we excluded the work of INisini et al.l ()2010l ). because 
they varied Tmin ioi the model application. 

Table [5] shows that n(H2) and b vary with the esti- 
mated levels {v = 0, J). ri(H2) and b are both tend to be 
smaller when the model applied to lower J levels; n(H2) 
shows a more drastic variation than b . This tendency was 
already noted in Shinn et aj] (|2010D through the model 
apphcation to the data of the SNR HB 21 and OMC- 
1; now, the tendency is strengthened by additional data 
from other SNRs and outflows of YSOs. This indicates 
that the level population of shocked H2 gas is not de- 
scribed by the power-law admixture model with single 
n(H2) and b. Moreover, it suggests we should analyze 
the level population over as many levels as we can get; 
if not, we may have a limited picture on the shocked H2 
gas. 

In the point of extending the data coverage, our obser- 
vations are important because we extended the v — 
level population obtained from previous mid-infrared 
studies, from E{v,J) < 7,000 K up to < 22,000 K 
(see Fig. H]). We could reproduce the estimated H2 
population with two n(H2) and one 6 after including 
the H atom as an additional collision partner (Table 
HI). The derived two n(H2)s, n(H2)=102-^~^ ® cm'^ and 
n(H2)=10^ '*^^'^ cm^"^, fall into the range previously ob- 
tained, n(H2)=10^ '^-'^ ° cm-^ (Table E]). As distinct 
from the previous applications (Table [S|), we included H 
atoms in the model, which efficiently excite high- J v — 
levels; hence, we obtained smaller n(H2) for the higher 
density component than we did without H atoms, by a 
factor of 0.8 dex and 1.7 dex for the clump C and G, 
respectively. 

As mentioned in section the lower and higher den- 
sity component mainly contribute to the lower and higher 
energy levels, respectively (Fig. ITO)) . This clearly show 
why we get different n(H2)s depending on which levels 
are used for the estimation (cf. Table [5]). This den- 
sity property also indicates that there may be a non- 
uniqueness in the model parameter, because the observed 
population constrains only small portion of the modeled 
population (cf. section 13.41) . Indeed, the T^ax of the 
lower density gas and the Tmin of the higher density gas 
can be lower and higher, respectively, since small varia- 
tions of those values may have a negligible effect on the 
fitting results. We checked these possibilities by vary- 
ing Tmax and Tmim with the obtained fitting parameters 
(Table |4|). In the case of clump C, we found that the 
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Tmax of the lower density gas cannot be decreased from 
the initial model setting 4000 K, while the Tmin of the 
higher density gas can be relaxed up to 1000 K. This 
result means that a small amount of warm, high den- 
sity H2 gas can explain the high energy-level population, 
i.e., iV(H2; T > 1000 K)~ 4 x 10^^ cm'^. In the similar 
sense, the Xh in the lower density gas can be lower, even 
down to the hydrogen-free case, n(H I)/n(H2)=0, since 
the excitation of the low-J v = levels is dominantly 
determined by the collisions with H2 rather than H I. 

4.1.2. Nature of the Shocks 

There have been many observational studies trying to 
identify the shock type at the clump C and G at diverse 
wavelengths — millimeter, sub-millimeter, and infrared. 
Much of the observational resul ts were compared with 
several J- and C-type shocks bv lSnell et al.l ()2005D . and 
they concluded that a combination of shocks (dissociative 
and non-dissociative) is required to explain all the ob- 
servational results. Later, Spitzer mid-infrared spectral 
observations showed the emission lines likely from such a 
combination of dissociative (J-type ) and non-dissociativ e 
(C-type) shocks at the clump C (|Neufeld et al.l I2007D . 
Besides, previous studies on the shocked H2 gas of pro- 
tostellar outflows which covered a similar E{v, J) range 
with ours (~ — 25, 000 K) showed that a mixture of 
C- and J-type p roperties are required to explain the H2 
level population (iFlower et al.l2003HGiannini et al.ll2006l: 
IGusdorf et al.ll2008n . [ 

Under the shock-combination preference, the immedi- 
ate question is which type of shocks the two component 
H2 gas originate from. We first considered the lower den- 
sity gas which occupies most mass of the shocked H2 gas 
(Fig. [TU] and Table 0]). If the lower density gas originates 
from J-type shocks, the recombinational lines of H I are 
also expected. We checked the relative intensity of Bra 
4.05 /im line a nd H2 = — S(3) line from a the- 
oretical model (jHollenbach fc McKed Il989| ). The ratio 
of Bra/[H2 v = - S(3)] is between ~ 1 - 10 over 
the shock velocity 50 — 100 km s~^ and preshock density 
10^ - 10^ cm-3.' 

The H2 density obtained from the model fitting likely 
represent the postshock density, since the II2 emissions 
are mainly e manating from the reforme d H2 gas in J- 
type shocks (|Hollenbach fc McKed 1 19891 ). The typical 
J-type shock velocity is around 100 km s~^ under a 
gener al interstellar medium condition (|Draine fc McKed 
Il993f l , hence the density of the molecular reformation re- 
gions at far down stream would be > 100 times higher 
than the preshock density; this density r atio can be 
seen from many numerical studies (e .g. lAllen et al.l 
[2OO8I: IFlower fc Pineau Pes Foretsll20"l0l ). Therefore, the 
preshock density should be n(H2) < 10^'^ cm~'^ for 
the lower-density gas to be from J-type shocks. This 
preshock density range is somewhat lower than consid- 
ered in the above theoretical model (10"^ — 10® cm~^), 
but the similar ratio BrQ!/[H2 v — — S(3)] would be 
maintained since the postshock structure is insensitive to 
the preshock density in J-shocks. Therefore, if the lower 
density gas is from J-shocks, then there should be Bra 
of ~ 2 X 10^'^ — 2 X 10~^ erg s~-^ cm~^ sr~^ based on the 
observed intensity of H2 u = — S(3) ()Neufeld et al.l 
l2007i) . However, no Bra line of such a high intensity was 



detected (Fig.©. 

In addition, theoretically, a J-shock produce more 

abundant high-T (~ 10'^ K) H2 gas than a C-shock 

JWilgenbus et al.l l2000t IFlower fc Pineau Pes ForetsI 
2010). If the lower density gas is from J-shocks, the 
H2 population diagram must show a much up-turn 
curvature than we obtained (Fig. irU]) . Overall, the lower 
density gas is not likely from J-type shocks. 

We then checked whether the lower density gas orig- 
inate from C-type shocks. The postshock gas behind a 
planar C-type shock can be approximated as an isother- 
mal gas whose density n(H2 ) and temperature Ts are as 
follows (|Neufeld et alJl2006li . 

n(H2) = 1.5no (3) 
T,^375bg°'"'vli^ K, (4) 

where 10us6 km s~^ is the shock velocity and 
^^(no/cm"^)^/^ /iG is the assumed preshock magnetic 
field for a shock propagating in material of preshock H2 
density uq. We here note the abo ve approximation wa s 
tested over no = 10"* - 10® cm'^ (|Neufeld et al.ll2006D . 
From this approximation, we can say the lower den- 
sity gas is from C-shocks propagating into a preshock 
medium ri(H2) ^ iq2-6-3.6 qj-^-s^ Interestingly, this den- 
sity is similar with the typical density of molecular clouds 
(> 10"^ cm""^), where th e hydrogen dominantly exists in 
the molecular form H2 (jSnow fc McCallll2006f ). Such a 
similarity gives a natural way to explain the mass domi- 
nance of the lower density gas (Fig. fTO)) . since H2 survives 
C-type shocks. Therefore, C-type shocks are more suit- 
able to explain the lower density gas than J-type shocks. 

If the lower density gas is from C-type shocks, the 
power-law index b may reflect the multiple shocks of dif- 
ferent velocities. According to the approximation with 
6b = 1, it ranges from 3 to 58 km s~^, which corre- 
sponds to ^ 100 — 4000 K gas. This velocity range is 
theoretically expected for a C-type shoc k whose preshock 
densi ty is ^ 10'^ cm~'^ and 6_b = 1 (|Le Bourlot et al.l 
|200^ . Also, such a diversit y in the shock veloc ity is 
consistent with the results of iHewitt et al.l ()2009aD who 
found two C-shocks are preferred to explain the exci- 
tation of shocked H2 gas around SNRs over the range 
E{v, J) ~ — 8, 000 K. However, we here note that this 
shock velocity range depends on the model parameters 
we fixed {Tmin, Tmax)] therefore, the shock velocity range 
mentioned above (3 — 58 km s~^) is not observationally 
determined one. Instead, it should be understood that 
multiple shocks of a few to a few tens of kilometer per 
second are required to explain the level population of the 
lower density gas. 

The diversity in the shock velocity may be a natural re- 
sult during the shock-cloud interaction, since the shock 
driving pressure would be higher at the head than the 
side. This kind of property can be inf erred from a nu- 
meric al simulation for radiative clouds (jNakamura et all 
|2006[) . which showed that the velocity dispersion of the 
shocked cloud is larger along the shock-normal (the blast 
wave direction; ^-axis iu iNakamura et al. 200 6) than the 
shock-tangential (ro-axis in lNakamura et aL[l2006t ): this 
velocity dispersion does not exactly corresponds to the 
dispersion of the shock velocity in the cloud, however, 
it shows the tendency which direction the shock would 
propagate faster and cause the larger velocity dispersion. 
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We assigned the same H I abundance both for the lower 
and higher density gas, and obtained Xh = —1-7, equiv- 
alent to n(H I)/nH = 0.01 and n(H I)/n(H2) = 0.02 (Ta- 
ble HI . This value is higher for typical molecular clouds 
(riH > lO"* cm-3, iDraine et all 11983!). However, as men- 
tioned in section 14.1.11 we can achieve the population 
of the lower density gas even with zero H I abundance, 
since the population is dominantly determined by the 
collisions with H2 rather than H I. Considering this fact, 
the C-shock interpretation for the lower density gas has 
no contradiction with the H I abundance. 

Now we move to discuss which type of shock the higher 
density gas originate from, considering the H I abundance 
first. The H I abundance, Xh^ is mainly determined 
by the vertical gap between u = and u > levels in 
the population diagram, and the gap is dominantly de- 
scribed by the higher density gas (see Fig. [TOl) . The ob- 
tained Xh = —1.7, equivalent to n(H Ij/rin = 0.01 and 
n(H I)/n(H2) — 0.02, is expected for a typical diffuse 
cloud (uh ~ 10^ cm~'^), but high for a typical molec- 
ular cloud (tt - h ^, 10^ cm~'^) which has Xh < —4.0 
(jPraine et al.|[T983h . Therefore, the obtained Xh must 
be from one of the following cases: (1) non-dissociative 
(C-type) shocks propagating into molecular clouds of 
the high Xh value; (2) non-dissociative (C-type) shocks 
propagating into diffuse clouds; (3) dissociative (J-type) 
shocks propagating into molecular clouds. 

We c hecked the first cas e. In view of the c alcula - 

tions mon fc Werned (fT97l . lDraine"etaI1 (fT983h . 

and [Goldsmith fc Lil (|2005[ ). the total ionization rate 
for H2 must be C ~ 10"^^ - 10"^^ s'^ to make the 
obtained abundance in the higher density H2 gas of 
n(H2) ~ 10^" cm~^. However, the measured rate 
is m uch smaller than the required: C ^ 5.5 x 10^^^ 
s-^ (IHewitt et al.l l2009bD and C ~ 2.6 x lO-^^ s"! 
(jlndriolo et al.ll2010[ ). The second case is negative either. 
It was shown that the postshock gas behind a planar C- 
type shock can be approximated with a i sothermal gas 
whos e density is equal to ~ l-5npreshock ()Neufeld et al.l 
I2OO6D . If the shock is C-type, the obtained density 



of n(H2) = 10 



5.4-5.8 



cm is contradict with the re- 



quirement that the preshock medium must be a dif- 
fuse cloud (tt-h ^ lO'^ cm"'^). In the third case, how- 
ever, the high Xh of —1.7 is expected under appropri- 
ate conditions; such an Xh can be achieved at the T ~ 
1000—4000 K region where t he H2 emission is mainly em- 
anating fr om reformations ("Hol lenbach fc McKed fl989l : 
iFlower fc Pineau Dcs Forcts 201^^ 

If the higher density gas is from J-type shocks where 
the H2 emission dominantly comes from the postshock 
reformation regions, the obtained model parameters for 
the higher density gas (Table S]) can be regarded to rep- 
resent characteristic physical parameters of the reforma- 
tion regions. Based on this, the preshock density should 
be a factor of < 100 smaller than the obtained density 
of n(H2) = lO^ "'^^-^ cm^'^ since the typical velocity of 
J-shocks is around 100 km s~^. Also, there should exist 
shocked H I gas, the raw material for the reformation of 
H2. The shocked H I gas has already been observed along 
the " W" ridge (Fig. [T| ) we are studying (Braun fc Stroml 
ll986HLee et"alll2008ari . Eei (|2007D measured the column 
densities of the shocked H I gas at the clump C and G, 
which are about ~ 3 x 10^^ cm^^. If we assume that the 



cm is obtained, which is smaller 

5.4-5.8 cj^-3_ 



shocked H I spreads over along our line-of-sight as long as 
the width of the "W" ridge (^^ 2' '--^ 1 pc), then the mean 
H I density of - 10 

than the postshock H2 density ri(H2) — 10 
These n(H I) and n(H2) are consistent with the J-shock 
interpretation, since the H I recombines before the H2 
and the postshock density increases along the distance 
from the shockfront. We also ch ecked the line profile 
of H I 21cm emission (lLea[200l and H2 u = 1 ^ 
S(l) emission (jRosado et al.l 120071 ) at the clump C, and 
found that they both show the shocked component at the 
similar velocity ~ —30 km s~^. Therefore, the J-shock 
interpretation for the higher density gas seems plausible 
again. 

We then considered the meaning of the power-law in- 
dex b in the higher density gas. As mentioned in section 
14.1.11 Tjnin of the higher density gas can be increased 
from 100 K to 1000 K without making any difference on 
the final level population of E{v, J) > K levels. This 
means that the gas in T = 1000 — 4000 K is necessary 
to describe the observed population of E{v, J) > lO'' K 
levels. We thus checked whether the b variation makes 
any significant change to the level population of E{v, J) 
> 10'' K, summed over T = 1000 - 4000 K; the obtained 
n(H2) and Xh was fixed while the b value was varied. 
We found that the b variation makes a negligible effect 
over b — 1.0 — 4.0. This is because the population of 
E{v, J) > lO'' K levels has an almost same slope over 
the temperature integration range (1000 — 4000 K) in 
the power-law thermal model. Therefore, the b value has 
almost no practical meaning for the higher density gas. 

From all the above discussion, the lower and higher 
density gases are likely from the C-type and J-type 
shocks, respectively. This conclusion is consistent with 
the conclusion from previous studies which claimed 
a mixture of dissociative and non-dissocia tive shocks 
for the description of the observed results (jSnell et al.l 
120051: iNeufeld et alll2007[) . If we interpret our results 
based on t he hierarchical pict ure of molecul ar clouds 
(|Bergin fc T afalla 2007.; Will iams et al.l [20Ciol ). the C- 
type shocks propagate into "clumps," while the J-type 
shocks propagate into "clouds" (interclump media); the 
typical densities of these two constituents ( "clouds" and 
"clumps") are 50 — 500 cm~^ and lO'^ — 10^ cm"'^, respec- 
tively. This interpretation is consistent with the require- 
ment that the ram pressure (^^ pv'^) should be similar 
between C-type and J-type shocks. In the above dis- 
cussion, we claimed that C-type shocks propagate into a 
preshock medium of ^ 10'^ cm~^ with a shock velocity 
of a few 10 km s~^, while J-type shocks propagate into 
a preshock medium of < 10'^ cm"'^ with a shock velocity 
of ~ 100 km s~^ 

Before closing this section, we note that the H2 w = 
1 S(l) is dominantly originating from J-type shocks 
rather than C-type shocks, according to our results. In- 
terestingly, it was observed that the [Si II] 34.8 /im emis- 
sion well foll ows the H2 v = 1 S(l) emission a t 
the clump C ()Richter et" al.lll995at INeufeld et al.ll2007D : 
the [Si II] 34.8 /um emission is an efficient cooling line 
at T < 5000 K where the H 2 reformation also happens 
(jHollenbach fc McKeeiri989f ). The intensity of [Si II] and 
H2 w = 1 S(l) are 0.97 and 1 8 in 10^"^ erg s"! cm"^ 
sr-\ respectively (jNeufeld et al.ll2007l Table EJ. The ra- 
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tio [Si II]/H2 w = 1 -> S(l) is 0.54, which is significantly 
lower than the theoretical expect ation of a 100 km 
J-shock into 10^ cm"^ gas, ~ 5.0 (|Hollenbach &: McKe^ 
[1989) . In spite of this difference, the J-shock origin of 
H2U = 1^0S(1) emission is worth to check further 
whether it is general in shocked regions, considering that 
the shocked H2 gas usually shows a similar shape in the 
level population over E{v, J)= — 25, 000 K (cf. section 
14X31) . 

4.1.3. Embedded Coherence in the Obtained Model 
Parameters 

Table |4] shows some coherences of the obtained model 
parameters for the two-component gases, although the 
clump G has no data for the population of E{v, J) — (0,2) 
and (0,3) which are important in the determination of b 
and density. For example, the density ratio and the col- 
umn density ratio of the two-component gases are not 
much different at the each clump, C and G. This may 
originate from some common proper ties of the pre shock 
medium. One step forward, Riclitcr et al.l (jl995ai rbf) al- 
ready noted that the similarity between the H2 level pop- 
ulation of SNRs and outflows of YSOs. 

We can study such similarities with the two-ri(H2) 
power-law thermal admixture model by comparing the 
model parameters, and may lead out the properties of 
the preshock medium. Especially, the power-law index b 
would give information about how shock velocity is di- 
verse at the shock-cloud interaction regions, which may 
be related with the geometrical structure of the clouds. 
The coming James Webb Space Telescope which can per- 
form spectral observations over the ~ 1 — 30 /im wave- 
length range would provide an excellent opportunity for 
such studies. 

4.2. Warm Diffuse Background H2 Gas 

At the background region (BG in Fig. [1]), we de- 
tected the H2 u = 0-0 S(9) hue (Fig. ^. The 
extinction corrected column density gives log[N(H2; 
u, J=0,ll)]=15.08 ± 0.11 cm-2 (Tabled. This emis- 
sion may be related with the diffuse background H2 gas , 
previously observed towards SNRs (jNeufeld et al.l [20071 : 
iHewitt et al.ll20d9al ). In order to check whether the ob- 
tained (u, J=0,11) population is plausible for such an 
diffuse background H2 gas, the population was compared 
with those obtained towards other Galactic lines of sight. 
Since the relative shape of the level population deter- 
mines the physical condition of the II2 gas, we must com- 
pare such an relative shape, thus we need the population 
of other levels in addition to (w, J=0,11) level. 

For this, we employed the (u, J=0,2) population of 
the diffuse background II2 gas, previously observed to- 
wards the clump C of IC 443 (Ncu fcld ct aL. 2007), as the 
(u, J=0,2) level population of our BG region (Fig. [1]). 
We adopted Idiff = 1-5 x 10~^ erg s~^ cm~^ sr^^ as 
th e intensity for t he H2 w - S(0) line (cf. Fig. 19 
of iNeufeld et al.l l2007f ). and corrected the extinction 
with Av = 13.5 e mploying the "Milky Way, -Rv=3.1" 
extinction curve (jWeingartner fc Draind 120011 : iDraind 
[200I . The population is log[N(H2; v, J=0,2)]=19.05 ± 
0.07 cm~^. After all, the weighted populations of the 
(u, J)=(0,2) and (0,11) levels for the BG region are 
log[N(i;, J)/gj]=18.3 and 13.2, respectively. 



We compared the populations of these two levels with 
those of (u, J)=(0, < 8) obtained toward the "translu- 
cent lines of sight" to the Galactic background stars, 
from the far-ultrav i olet o bservations of II2 absorption 
lines (jJensen et al.l 12010( 1. Since the p opulation of 
V, J) =(0, 11) is absent in the results of 'J ensen et al.l 
2010[) , the relative curvature in the population diagram 
was compared, by scaling up or down the populations of 
(u, J)=(0,2) and (0,11) levels. From the comparison, we 
found that our populations make a more upturn curva- 
ture than expected from the excitation tempera ture of 
J > 2 levels, ~ 200 - 550 K (jJensen et al.ll2010D . This 
means the existence of H2 gas with higher temperature 
than 550 K, which is consistent with the suggestion of 
I Jensen et an()2010[ ): more than two temperature compo- 
n ents are probably ne cessary to describe the population. 

iJensen et al.l (12010) tried to fit the population obtained 
toward two sample targets, with a model for photodis- 
sociation regions; however, the fi ts were poor. Su ch a 
mismatch was a l so rec ognized by iGrv et all ()2002D and 
iFalgarone et al.l ()2005[ ). and they suggested that another 
collisional excitation mechanisms are required to explain 
the observed population, such as magnetohydrodynamic 
shocks or intermittent dissipation of turbulence. The 
diffuse background II2 gas towards the SNR IC 443 may 
originate from these mechanisms. 

5. CONCLUSIONS 

We present near-infrared (2.5 — 5.0 /im) spectra toward 
three shocked H2 clumps (B, C, and G) of the SNR IC 443 
and one background (BG) region (cf. Fig.[T]). The obser- 
vations were performed with the satellite AKARI during 
its post-Helium phase. Only II2 emission lines were de- 
tected toward all four directions. We measured the line 
intensities and, after the reddening-correction, obtained 
the relevant level populations. The level populations 
were compared with the ones from previous near- and 
mid-infrared observations. AKARI level populations are 
well fitted with those from previous mid-infrared space 
observations, while there is a systematic difference with 
those from previous near-infrared ground observations. 
We attributed this difference to the calibration error of 
previous near-infrared observations, although it may be 
caused by the different position of observed regions. 

With the AKARI near-infrared observations, we could 
extend the level population of shocked H2 gas in the 
clump C and G obtained from mid-infrared observations, 
from 7000 K to 22000 K continuously. From these 
combined level populations, we found that the v = lev- 
els are more populated than the v > levels at a fixed 
level energy, which means the population cannot be re- 
produced with any combination of H2 gas in LTE, usually 
used for the description of shocked II2 gas. The popula- 
tions were well described two-n(Il2) power-law thermal 
admixture model, including the H atom as a collisional 
partner in addition to H2 and He. The model parame- 
ters are two number densities n(H2), two column densi- 
ties A^(H2), one power-law index b, and the H I to H2 
abundance Xh (cf. Table [4]). 

We attributed the lower (n(H2) ~ xo^-^~^-^ cm""^) and 
higher (n(H2) = 10^ '*~^-* cm^'^) density gases to the 
shocked H2 gas behind C-type and J-type shocks, respec- 
tively, based on several arguments such as the preshock 
density, the amount of shocked gas, the obtained high 
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H I abundance, the shape of level populations, the line 
profile, etc. This interpretation is consistent with the hi- 
erarchical picture of molecular clouds whose constituents 
are classified into "cl o uds," "clumps , " an d "cores" 
(|Bergin fc Tafallal [2007t IWilliams et ai1[20Q0l) . The C- 
type shocks are likely propagating into the "clumps" 
10'^ cm~^), while the J-type shocks are propagat- 
ing into the "clouds" (interclump media, < 10^ cm~^). 
The power-law index b, mainly determined by the lower 
density gas, is attributed to the diversity in the shock 
velocity propagating into the clouds. Such a velocity di- 
versity may be a natural result during the shock-cloud 
interaction, since the shock driving pressure would be 
higher at the head than the side. According to our re- 
sults, the H2 f = 1 — > S(l) emission is mainly from 
J-type shocks propagating into interclump media. In ad- 
dition, our power-law thermal admixture model would 



be useful to perform statistical studies on the shocked 
II2 gas observed around other SNRs and YSOs, which 
shows coherent excitation characteristics. 

H2 emission was also detected at the BG region, and 
we attributed it to the diffuse H2 gas, pervaded in the 
Milky Way. This gas may be excited by collisional pro- 
cesses, like shocks or turbulence dissipation, in addition 
to ultraviolet photon pumping. 

This work is based on observations with AKARI, a 
JAXA project with the participation of ESA. The au- 
thors thank all the members of the AKARI project. J. 
H. S. is grateful to Jae-Joon Lee for the useful discus- 
sion on the two-density interpretation. This research has 
made use of SAOImage DS 9, developed by Smith sonian 
Astrophysical Observatory (jJove fc Mandell [20031 ). 

Facilities: Akari 
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Fig. 1. — The observed slit positions. Slit positions are indicated on the 2MASS Ka image of the SNR IC 443 HSkrutskie et al.ll200g '). Slit 
positions are named after their representative clump names — 'B', 'C, and 'G.' 'BG' means the background. The red 'W ridge is dominated 
by H2 emissions. 
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Fig. 2. — The AKARI IRC near-infrared spectra. Singular and blended lines are indicated with vertical bars and crosses, respectively. 
These lines are all the H2 emission lines (cf. Table|2]|. In the case of 'BG', a twenty-times enlarged spectrum is plotted for the identification 
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Fig. 3. — Fitting for the H2 emission lines, observed toward the clump B. The solid and dashed lines indicate the continuum+line and 
continuum, respectively. The blended line components are indicated by the dotted lines, if any. The feature around 2.56 /im is ignored 
since it is near the end of the band coverage. 
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Fig. 4.— 



Fitting for the H2 emission lines, observed toward the clump C. The rest is the same as Figure [3] 
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Fig. 4. — Continued. 
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Fig. 5.— 



Fitting for the H2 emission lines, observed toward the clump G. The rest is the same as Figure [3l 
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Fig. 9. — Continued. 
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Fig. 10. — The two-component model fitting results for the H2 population observed in the clump C (top panels) and G (bottom panels). The 
line-connected gray points in the left panels are the total H2 population obtained from the model fitting, and their individual components 
are displayed in the right panels as red (lower density) and blue (higher density) ones, (see text for the model description.) The rotational 
quantum number (J) is printed out near the corresponding point. 
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Fig. 11. — The contour of in the plane of model parameters for the clump C (left panels) and G (right panels). The 68% and 95% 
confidence levels are outlined. The tick marks along the contours indicate the directions that values are decreasing. The 'ffi' indicates 
those model parameter values, whose is minimum, i.e. the best fit. In the right-bottom panel, the x^ values are plotted about the 
parameter b only, since the Xjj is fixed (cf. Table|4l the bottom figures are displayed in the following page.). The dotted-line indicates the 
best-fit parameter, and the dashed-lines indicate the 68% and 95% confidence levels. 
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TABLE 1 

Summary of the AKARI IRC Observations 



Region 



Pointing Position 
(RA, Dec; J2000) 



Observation ID AOT'' 



B (06:17:16.3, +22:25:41.0) 1420803-00 

C (06:17:44.2, +22:21:49.1) 1420804-00 

G (06:16:41.8, +22:31:41.0) 1420805-00 

BG (06:17:30.0, +22:17:00.0) 1420806-00 



IRCZ4 
IRCZ4 
IRCZ4 
IRCZ4 



" Astronomical Observation Template. It is a pre-defined obser- 
vation se quence. See th e IRC Data User Manual for Post-Helium 
Mission. jOnaka et al.ir200a) . 



TABLE 2 

Observed H2 Emission Lines toward Each Region 



Transition Wavelength Observed Intensity 

(/im) (10~® erg s^^ cm^^ sr^^) 













B 


C 


G 


BG 


V 




1-0 


0(2)=^ 


2.63 


<37.9±6.1 


<93.4+10.4 


<58.2+7.1 




V 




2-1 


Q(9)'' 


2.72 




<25. 1+5.3 


<13.5+4.2 




V 




1-0 


0(3) 


2.80 


89.9+10. 


323.1+32. 


221.9+22. 




V 




2-1 


0(3) 


2.97 


<13.5 


27.0+4.8 


18.5+4.1 




V 




1-0 


0{4) 


3.00 


39.8+5.2 


96.7+10. 


63.0+6.9 




V 




1-0 


0(5) 
0(6)^ 


3.23 


61.1+7.4 


196.2+20. 


127.0+13. 




V 




1-0 


3.50 


<18. 2+4.1 


<46.5+5.8 


<29. 2+4.3 




V 




0-0 


S(15) 
S(14)'' 


3.63 




47.8+5.6 


24.1+3.4 




V 




0-0 


3.72 




<26.8+3.9 


<10. 7+2.6 




V 




1-0 


0(7) 


3.81 


18.9+3.4 


54.1+6.1 


36.7+4.2 




V 




0-0 


S(13) 


3.85 


27.6+3.9 


96.0+10. 


53.1+5.7 




V 




0-0 


S(12) 


4.00 




38.5+4.6 


21.5+3.1 




V 




1-1 


S(13)'' 


4.07 




<24.1+3.5 


<14.9+2.7 




V 




0-0 


S(ll) 


4.18 


60.0+6.8 


199.3+20. 


113.4+11. 




V 




0-0 


S(10) 


4.41 


40.1+5.2 


122.8+12. 


71.2+7.6 




V 




1-0 


0(9) 


4.58 




16.0+3.2 


9.6+2.6 




V 




0-0 


S(9) 


4.69 


134.4+14. 


430.8+43. 


265.7+27. 


14.2+3.5 



Note. — For those lines whose significance is lower than 3.0, the intensities are ex- 
pressed with 90% upper confidence limits. The error includes both statistical and sys- 
tematic components. See text for detail. 

^ These lines may be blended with nearby lines, hence we indicated with the '<' sign. 
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TABLE 3 

Extinction Corrected H2 Column Density toward Each Region 



State 


Energy Level 




log N(H2 


; ^- J) 




{v,J) 


(K) 




(cm~ 










B 


C 


c 


BC 


f r\ 1 1 \ 
yyj^LL ) 




ID.UDitU.Uo 


iD.0DitU.U4 


iD.ozitu.U4 


1 c: nQ_l_n 1 1 
iO.UoitU.ii 


(f\ 1 o^ 


1 iy'-tu. 


iO.oDltU.UD 


io.oDitU.U4 


io.ooitU.UO 




( f\ 1 '^^ 
^U,io j 


lo / Uo. 


i0.4UltU.U0 


io.y^itU.U4 


io.D4itu.U4 




(c\ 1 A\ 


1 c:c: /in 




io.uomu.uo 


1 4 7Q_i_n c\P, 




(f\ 1 K\ 

l^U,iD ) 


1 7/1 /I Q 


1 /I QQ_i_n nfi 


10. / itU.UD 


10. U ( itU.UO 




(0,16)^ 


19403. 




<14.72±0.06 


<14.28±0.10 




(0,17) 
(1, 0)=^ 


21411. 




14.89±0.05 


14.55±0.06 




5987. 


<15.31±0.07 


<15.70±0.05 


<15.41±0.05 




(1, 1) 


6149. 


15.96±0.05 


16.52±0.04 


16.28±0.04 




(1> 2) 


6471. 


15.76±0.06 


16.14±0.05 


15.89±0.05 




(1. 3) 


6951. 


16.07±0.05 


16.58±0.04 


16.33±0.05 




(1, 4)=^ 


7584. 


<15.68±0.10 


<16.09±0.05 


<15.84±0.06 




(1> 5) 


8365. 


15.85±0.08 


16.30±0.05 


16.09±0.05 




(1, 7) 


10341. 




16.12±0.09 


15.87±0.12 




(1,15)=' 


22516. 




<14.84±0.06 


<14.60±0.08 




(2, 1) 


11789. 


< 14.94 


15.25±0.08 


15.01±0.10 




(2, 9)=^ 


18107. 




<15.56±0.09 


<15.21±0.14 





Note. — The extinctions arc corrected, employing the "Milky Way" extinction curve 
{Rv — 3.1. [Weingartn er &: DrainQii2001i), with these parame ters: Av — 13.5 for clump B, 
clump C and BG |Burton et al.[|1988l : l>Jeiifeld &: Yuaj3l2008D and Av ^ 10.8 for clump G 
ijRichter et al.lll995al ). respectively. 

^ The population cannot be not determined because of the probable line blending with 
nearby hnes. cf. Table [2] 



TABLE 4 

Fitting Results for the Model Parameters 





component 1 


component 2 








Region 


Iog[Af(H2)J 


logln(H2)J 


Iog[Af(H2)| 


log[n(H2)| 


b 








(cm~^) 


(cm"-^) 


(cm~^) 


(cm~3) 






(= xVd.o.f) 


IC443C 
IC443G 


21.0^0;? 

OQ r)"I"t).2 

^^•^-0.4 


9 o+O.l 

o O+0.2 
3.»_0.4 


IQ 4+0.3 
21 l+O-** 


5 4+0.3 

cr O + 0.8 
O.S_0.4 


1 6+°-3 
o tr+0.3 
'^•'^-0.6 


1 7+0.2 
'-0.2 

-1.7=" 


5.1 (=66.1/13.0) 
3.4 (=40.6/12.0) 



Not e. — The confidence limits are given with a 95% significance (cf. Fig. IIH . A^(H2) means iV(H2; T > 100 K). See 
section |3.4| for the detailed description about the parameters. A'^(H2) of the component 2 is not firmly determined value, 
since it sensitively depends on the Tj^in of the model (see section |4.1.H . 

^ In this case, the parameter Xh is fixed as —1.7, to increase the degree of freedom. See SGction |3.4| for more detail. 



TABLE 5 

Obtained n(H2) and b from Previous Studies 



Target 


log[n(H2)] 
(cm ■^) 


b 


Instrument 


Estimated Levels 


Ref. 


SNR IC 443 


5.0 - 7.0 


3.0 - 6.0 


Spitzer IRAC 


1, = 0, J = 6 - 


15 


1 


SNR HB 21 - Cloud N 


2.7-3.3 


2.9 


AKARI IRC 


u = 0, J = 3 - 


8 


2 


SNR HB 21 - Cloud S 


4.6 


4.2 


AKARI IRC 


v = 0,J = A- 


13 


3 


OMC-1 


3.5 


1.9 


AKARI IRC 


v=0,J =3- 


8 


3 


OMC-1 


6.0 


3.7 


AKARI IRC 


i, = 0, J = 4- 


13 


3 


outflows of YSOs 


3.5 - 3.8 


2.3 - 3.3 


Spitzer IRS 


^; = 0, J = 2 - 


9 


4 


outflows of YSOs 


^ 7.0 


3.0 - 5.5 


Spitzer IRAC 


u = 0, J = 6 - 


15 


5=^ 


outflows of YSOs 


^ 5.0 


3.0 - 6.0 


Spitzer IRAC 


u = 0, J = 6 - 


15 


6 


SNRs and outflows of YSOs 


3.3 - 3.6 


2.3 - 3.1 


Spitzer IRS 


v = 0,J = 2- 


9 


7 



Not e. — 1: Ncufcld & Yuan 12008), 2: Shinn ot al. (2009), 3: Shinn ct al. |20T3), 4: INeufeld et"all 
I I2009I) . 5: iLec ct al. (2010), 6: Takami ct al. (2010), 7: Yuan & Ncufcld |20lJ). The modeled gas 
consists of H2 and He only, with n(H2)=0.2 n{He). 

^ The central sources (position C and E) are excluded since H2 may not be the only contribution to the 
IRAC bands. 



